Observation of the T(2 3 Di) and indication of the T(l 3 Di) 
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We present evidence, from published BABAR data, for the existence of further excitations in the 
bb spectrum, namely the long-awaited T(l 3 Di) and T(2 3 _Di), with central masses of 10098±5 MeV 
and 10495±5 MeV, respectively The significance of the T(l 3 Di) and T(2 3 Di) signals is found to 
be 3.0 and 10.7 standard deviations, respectively. 
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In Ref. [l| the BABAR Collaboration presented 
an analysis of data on e + e~ — > 7r + 7r~ Y(25, 15) —> 
TT + n~e + e~ (and 7r + 7r~/i + /i - ), with the aim to study 
hadronic transitions between bb excitations and the 
T(15) and T(2S), based on 347.5 ftr 1 of data taken 
with the BABAR detector at the PEP-II storage 
rings. BABAR reported Q the first observation of the 
T(10580) — > r)T(lS) decay, with a branching fraction of 
S(T(10580) -> r?T(15)) = (1.96 ± 0.06 sta t ± 0.09 syst ) x 
10~ 4 , measured the value 2.41 ± 0.40 sta t ± 0.12 syst for 
the ratio of partial widths T [T(10580) -> t?T(15)] / 
r[T(10580) -> 7T7tT(15)], set 90%-confidence-level up- 
per limits on the ratios T [T(25) -> r?T(15)] / T [T(2S) 
-> 7T7rT(lS*)] < 5.2 x 10~ 3 and T [T(3S) -> ??T(15)] 
/ T[T(35) -> 7T7rT(15)] < 1.9 x 10" 2 and finally 
presented new measurements of the ratios T [T(10580) 
-> 7T7rT(2^)] / r[T(10580) ttttT{1S)} = 1.16 ± 
0.16 sta t ±0.14 syst and T [T(35) -> 7T7rT(2S)] / T [T(35) 
-»■ 7T7rT(15)] = 0.577 ±0.026 st at±0.060 syst . 

In the present paper the data of Ref. [l| are further 
analyzed for the possible existence of the T(l 3 Z?i) and 
T(2 3 Di) states, which we expect at roughly 2 x 4.724 + 
0.19 x (2n + 2 + 1.5) = 10.113 and 10.493 GeV @, for 
n = and n = 1, respectively. 

BABAR published Q the invariant-mass distributions 
of 34513 events for the reactions e + e~ — > n + iT~T(2S, IS) 
— > tt + ir~ fi + yT , and those of 53627 events for the reac- 
tions e + e~ — > Tv + ir~Y(2S, IS) — > ir + ir~ e + e~ . The events 
were preselected for the transitions e + e~ — > 7r + 7T _ T(lS') 
and e + e _ — > tt + t:~T(2S), according to the invariant 
masses of the final-state e + e~ and pairs. In the 

analysis, BABAR chose a rather large energy window 
around the T(15 l ) and T(25) masses, of the order of ±250 
MeV, for the invariant mass of the final-state e + e~ and 
pairs, but imposed further cuts on the sample of 
published data. In the present work, we shall apply much 
smaller windows, viz. ±20-30 MeV, though without im- 
posing more cuts. 

For a thorough analysis of the data we shall restrict 



ourselves to events from the BABAR sample of e + e~ — > 
tt + it~T(1S) — > 7r + 7r~e + e~ data. In Fig. Q] we show the 
invariant-mass distribution of the sample, for a window 
size of ±A (with A = 30 MeV) around the T(15) mass, 
and for a data binning of 10 MeV. Besides the dominant 
peaks of the T(25) and T(35), we observe a small en- 
hancement exactly where the data trend has a minimum 
due to the large signals. The enhancement is well visible 
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FIG. 1. Invariant-mass distribution for 2809 events (501 in 
the inset) in e + e~ — ► 7r + 7r~T(lS') —5- -K + Tv~e + e~. The data 
are taken from Ref. [l[. A window of A = 30 MeV around the 
T(15) is taken for the allowed invariant masses of the e + e~ 
pair. The bin width is 10 MeV. We have indicated our fit to 
the background by a solid curve. The T(l fli) signal has a 
relevance of 3. Oct. 

in the inset of Fig. [TJ 

In Fig. [1] we have also indicated, by error bars, the 
standard deviation on the number of events in each bin. 
For the data bins in the inset of Fig. (TJ we shall describe 
in detail how these values are determined. To that end, 
we study how the number of events in each bin grows 
when the window (A) is increased from to 200 MeV, in 
steps of 1 MeV. The result is depicted in Fig. [21 
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In order to clarify the relation between each set of data 
in Fig. [2] and the corresponding distributions shown in 
Fig. [TJ we have numbered the bins in the inset of Fig. [TJ 
as well as the sets of data in Fig. [2] Of course, each set of 
data in the latter figure should start out at events for 
A = (window closed). However, plotting the data this 
way would result in a confusion of points and curves. For 
clarity, we have thus increased the number of events in 
Fig. [2] by a constant, for each set of data. Hence, for all 
data sets the number of events in reality vanishes where 
the distributions meet the left-hand vertical axis. So the 
annotations on the axis denote the scale, and not the 
number of events corresponding to each data point. 
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we determine the size of the uncertainty in the number of 
events in each bin, also for the bins not shown in Fig. [2] 
One may easily verify that the sizes of the error bars 
shown in Fig. [JJ are in agreement with the deviations 
from smoothness in Fig. [3J 

One observes in Fig. [5] that the number of events in 
bin no. 7 increases faster for small values of A than the 
number of events in the neighboring bins. We believe 
this can be explained by assuming that events from the 
narrow-width T(l 3 Z?i) fall well within the 10 MeV bin for 
not too large windows, but start spilling to neighboring 
bins when the window is further opened. In the following, 
this issue will be dealt with more quantitatively. 

In Fig. [3]wc study, for the chosen signal bin (number 7 
in the inset of Fig. [TJ), the amount of events and also the 
signal-to-background ratio, as a function of the window 
size A. For that purpose, we determine the background 
signal by choosing the maximum of the number of events 
in the neighboring bins (numbers 6 and 8 in the inset of 
Fig. [J). 
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FIG. 3. Signal height (gray dots, pink in the online version) 
and signal-to-background ratio (black) for increasing window 
openings A from 1 MeV to 140 Mev in steps of 1 MeV, for 
bin no. 7 of the inset of Fig. [TJ Annotations on the vertical 
axis correspond to the signal-to-background ratio; the signal 
heights are 4 times larger. 



FIG. 2. Number of events for the reaction e + e~ — > 
7r + 7r _ T(lS) — ► tv + tt~ e + e~ , scaled up by a constant, for in- 
creasing window openings A, from 1 MeV to 200 Mev in steps 
of 1 MeV, in each of the 20 bins of the inset of Fig.fJJ The cor- 
responding bin numbers are given on the right-hand vertical 
axis. Each data set corresponding to a single bin is fitted by 
a curve quadratic in A (solid lines) that vanishes for A = 0. 
The data are taken from Ref. [l[. 

Each of the data sets corresponding to a particular 
interval in invariant e + e~ mass is fitted with a quadratic 
curve, where we assume that an infinite data set would 
lead to a smoothly growing number of events according 
as the window is further opened. From the standard 
deviation of the actual data with respect to smoothness, 



We have also studied all other bins shown in the inset 
of Fig. [TJ None of these has a behavior for the signal 
height and the signal-to-background ratio that is even 
comparable to what we observe in Fig.[3]for bin no. 7. For 
example, bin no. 12, at 10.15 GeV, has a maximum signal 
height of 8 events, which vanishes for windows larger than 
A = 32, while bin no. 13, at 10.16 GeV, has no positive 
signal for windows smaller than 80 MeV and, moreover, 
a totally negligible signal-to-background ratio for even 
wider windows. From all this we may conclude that, 
in the invariant-mass interval from 10.05 to 10.22 GeV, 
only one bin shows a very stable signal for all possible 
windows, and a reasonable signal-to-background ratio for 
windows up to about 70 MeV. 
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We assume that, for wider windows, possible events as- 
sociated with the supposed T(l 3 Z?i) signal start to spill 
over to the neighboring bins. As a consequence, the back- 
ground increases, whereas signal is lost. Hence, we prefer 
to study the statistical relevance for a smaller window, 
but still large enough for the signal height to be already 
reasonable. One may conclude from Fig. [3] that taking 
windows ranging from A = 25 MeV to A = 35 MeV is a 
good choice, since for lower values the signal is relatively 
small, whereas for higher ones the signal-to-background 
ratio decreases rapidly. For the purpose of our data anal- 
ysis, we have chosen A = 30 MeV, as shown in Fig. [T] 

The relevant bin for the T(l 3 £q) contains 30 ± 3.07 
events. For this bin, a global fit to the data (see 
Fig. [T]) gives a background of 18.0±2.53 events, leaving 
12.0±3.98 events for the signal. Thus, the signal has a 
significance of 3. Oct, which implies a strong indication of 
the T(l 3 L>i) at 10098±5 MeV. 

For the uncertainty in the mass of the T(l 3 £>i), we 
take half the bin width, as we do not expect this state to 
have a width larger than 1 MeV. Consequently, the scat- 
tering of the data within the 10-MeV-wide bin entirely 
stems from the method to collect data. 

The difference between the aforementioned value of 
10.113 GcV and the observed mass can be understood 
from the predicted small mass shifts of 3 Di bottomonium 
states due to open-bottom meson-loops [3L 

The CLEO Collaboration announced 4| the discov- 
ery of a J PC = 2~~ bb resonance with a central mass 
of (10161. 1±0.6±1.6) MeV. Although CLEO reported a 
significance of 10.2 standard deviations, their observa- 
tion is still today omitted from the Summary Tables of 
the Particle Data Group 0. Our identification of the 
1 — T(l 3 Di) state, which comes out 63 MeV below the 
2 T(1 3 D2), makes CLEO's assignment of quantum 
numbers very plausible, also in view of the quite differ- 
ent couplings of vector and tensor states to open-bottom 
channels. 

The observation of CLEO was recently confirmed by 
BABAR Q, obtaining a value of 10164.5±0.8 sta t±0.5 syst 
MeV for the T(l 3 I?2) mass. However, in the same paper 
BABAR claimed a mass of 10151 MeV for the T(l 3 Lq), 
which is close to the model prediction in Ref. [7[, but 
about 53 MeV above our value. In the present study, as 
mentioned before, we do not find any other relevant signal 
in the invariant-mass interval from 10.05 to 10.22 GeV. 

For the signal stemming from the T(2 3 Z?i), we follow 
precisely the same strategy as discussed above for the 
T(l 3 Di). In Fig.|4]we display how the number of events 
develops when the window (A) is increased from to 200 
MeV, in steps of 1 MeV, for all 10-MeV-wide bins shown 
in Fig. [6] The sets of data in Fig. |4l from the bottom 
to the top curve, correspond to the bins in Fig. [6] from 
left to right. For example, data set no. 12 counted from 
below in Fig. 21 which starts out at the annotation for 
500 events on the left-hand vertical axis, corresponds to 



bin no. 12 counted from the left in Fig. [6] This is the 
bin that contains more events than any other shown in 
Fig. [6] We observe for the latter set of data a substan- 
tially faster rise in the number of events for increasing 
values of the window size A, up to about A = 120, than 
for the data sets associated with the neighboring bins. 
The most obvious explanation for this behavior is the 
presence of an enhancement in bin no. 12, for which the 
most likely candidate is the T(2 3 Z)i). 
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of events for the reaction e + e — > 
7v + tt~ e + e~ , scaled up by a constant, for in- 



4. Number 
7r+7r-T(lS) 

creasing window openings A, from 1 MeV to 200 Mev in steps 
of 1 MeV, in each of the 34 bins of Fig. [6] Each data set cor- 
responding to a single bin is fitted by a curve quadratic in A 
(solid lines) that vanishes for A = 0. The data are taken from 
Ref. [J. 

In Fig. [5] we depict the signal and the signal-to- 
background ratio for the relevant bin. From this result 
we conclude that a window of A = 22 MeV is appro- 
priate for analysis. However, we admit that an equally 
good analysis is certainly possible for a window of A = 82 
MeV. We show the data for the latter window in Fig. [3 
using a data binning of 5 MeV, but shall not further anal- 
yse this case, for the same reason as we have discussed 
for the T(l 3 Di). Namely, for larger windows we suspect 
data to spill over over to neighboring bins. 
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FIG. 5. Signal height (gray dots, pink in the online version) 
and signal-to-background ratio (black) for increasing window 
openings A, from 1 MeV to 200 Mev in steps of 1 MeV, for 
the bin containing the Y(2 3 Di) signal of Fig. [6] Annotations 
on the vertical axis correspond to the signal-to-background 
ratio; the signal heights are 100 times larger. 



The data for the T(2 3 L>i) are depicted in Fig. El The 
relevant bin contains 130±4.51 events, while from a fit 
to the global data (solid curve in Fig. [S]) we find a back- 
ground signal of 81±0.81 events. For the signal we thus 
obtain 49±4.58 events, and so a significance of 10.7a. 
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FIG. 6. Invariant-mass distribution for 2311 events in e e~ — ► 
7r + 7r _ T(lS) — > ir + Tv~e + e~ , for window and bin sizes of A = 
22 and 10 MeV, respectively. Data are taken from Ref. 
Our fit to the background is shown by a solid curve. The 
T(2 3 Di) signal has a significance of 10.7<r. For completeness, 
we have also indicated the BB and BB* + BB* thresholds. 

Apparently, the 3 D\ bb states couple quite weakly to 
the transitions under study here, as can been seen from 
the large 2 3 Si and 3 3 Si peaks in Fig. [TJ to be contrasted 
with the more modest 2 3 Di (Figs. I6I7[) and especially 
1 3 Di (Fig. [IJ structures. This observation is in line 
with expectations from theory, since £>-wave quarkonium 
states are suppressed in e + e~ annihilation, due to the 
small wave function at short distances. However, the 



much stronger signal for the 2 3 Di state as compared to 
the 1 3 D\ seems puzzling at first sight, also in view of 
the 2 3 Si peak, which is higher than the 3 3 5i. A possi- 
ble explanation is offered by coupled channels. Namely, 
even bb systems below the lowest open-bottom threshold 
{BB) couple to virtual meson pairs, both the 3 S\ and 
3 Di states. This coupling inevitably gives rise to 3 5i- 3 L>i 
mixing. Now, the farther the bb system lies below the 
BB threshold, the more virtual will the meson pairs be, 
and so the smaller the 3 S\- 3 Di mixing. Therefore, the 
T(2 3 £>i), which lies only 64 MeV below the BB thresh- 
old, will have a significantly larger 3 S\ admixture than 
the almost 400 MeV lighter T(l 3 -Di), and so its effective 
coupling to e + e~ will be considerably larger. 

Nevertheless, the smaller 3 D\ signals may explain why 
they have not been discovered over the past three decades 
of e + e~ accelerator physics. Moreover, these states are 
probably also quite narrow. This complicates event se- 
lection, since on the one hand we prefer a narrow bin size 
for high resolution, but on the other hand we would also 
like to have sufficient statistics, which can be achieved 
by choosing large windows for the invariant masses of 
the final-state e + e~ pair around the T(15) mass. How- 
ever, too large windows inevitably produce several non- 
candidates and a very noisy background, as may be ob- 
served by comparing the distributions shown in Figs. [6] 
and [7] The choices we have made in the present analysis 
reflect a compromise between the two criteria. 
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FIG. 7. Invariant-mass distribution for 8337 events in e + e~ — ► 
7r + 7r~T(15') — ► 7r + 7r~e + e - for window and bin sizes of A = 82 
and 5 MeV, respectively. Data are taken from Ref. [J. For 
completeness, we have also indicated the BB and BB* + BB* 
thresholds. 

The here reported T(l 3 Di) and T(2 3 Di) masses lend 
further support to the level scheme of the harmonic- 
oscillator (HO) approximation to the resonance-spectrum 
expansion (HORSE) (see Ref. [§] and references therein) . 
In particular, the observed spacing of about 380 MeV be- 
tween the D states 0] makes a 3 3 D\ assignment for the 
T(10860) more plausible than 5 3 S 1 . 
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To be more precise, the pure HO, without accounting 
for meson loops, gives a mass of 2 x 4.724 + 0.19 x (4 + 2 + 
1.5) = 10.873 GeV 0, BABAR obtained 10.876+0.002 
GeV 0, while in Ref. @ we found 10.867 GeV for the 
T(10860) mass. Moreover, this result supports the con- 
clusion of our analysis in Ref. Q that the T(10580) is 
not the T(4ff) resonance, but rather a BB threshold 
enhancement [ll]. For the T(45) resonance, which in 
the HORSE is expected some 100-200 MeV below the 
T(3D), due to the effect of the meson loops, we found 
in Ref. [1] a central mass of 10.735 GeV, by analyz- 
ing BABAR data published in Ref. [13]. Finally, the 
conspicuity of the 2 3 Di signal, when compared to the 
much feebler 1 3 Di , is a compelling indication of virtual 
coupled-channel effects, as explained above. This repre- 
sents an additional endorsement of the HORSE and its 
nonperturbative description of decay, including open as 
well as closed meson-meson channels. 

In conclusion, we have observed excellent candidates 
for the T(l 3 Di) and T(2 3 Di) states, with central 
masses of 10098+5 MeV and 10495+5 MeV, respectively. 
Their observation, if confirmed independently, strongly 
supports the resonance level scheme of the harmonic- 
oscillator approximation to the resonance-spectrum ex- 
pansion for quarkonia. 
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